Diffractive deep inelastic scattering from multiple soft gluon exchange in QCD 
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Diffractive hard scattering is interpreted as the effect of soft gluon exchanges between the emerging 
energetic quarks and the nucleon's color field, resulting in an overall color singlet exchange. Summing 
multiple gluon exchanges to all orders leads to exponentiation and an amplitude in analytic form. 
Numerical evaluation reproduces the precise HERA data and gives new insights on the density of 
gluons in the proton. 
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Diffractive deep inelastic scattering (DDIS) in lepton- 
proton collisions involves hard scattering events where, 
in spite of the large momentum transfer from the 
electron, the proton emerges essentially unscathed with 
only a very small transverse momentum, keeping almost 
all of its original longitudinal beam momentum. The 
leading proton is well separated in momentum space, or 
rapidity y = 1/2 ln(£' +pz)/{E — Pz), from the central 
hadronic system produced from the exchanged virtual 
photon's interaction with the proton. Thus, this new 
class of events is characterized by a large rapidity gap 
(LRG) void of final state particles. 

Diffractive deep inelastic scattering (DIS) was discov- 
ered by the ZEUS and HI experiments at HERA 
but the first discovery of such a hard diffraction process 
was in pp collisions by the UA8 experiment These 
processes had actually been predicted Q by combining 
Regge phenomenology for low-momentum transfer (soft) 
processes in strong interactions via pomeron exchange, 
with large-momentum transfer (hard) processes based 
on perturbative QCD. By parametrizing the parton con- 
tent of an exchanged pomeron (or alternatively diffrac- 
tive parton density functions) it is possible to describe 
the HERA data. However, the extracted parton densities 
are not universal, since when used to calculate diffractive 
hard scattering processes in pp collisions at the Tevatron 
one obtains cross sections an order of magnitude larger 
than observed. 

As an alternative dynamical interpretation the Soft 
Color Interaction (SCI) model was developed in Ref. 
based on the assumption that the hard perturbative part 
of the interaction is the same as in ordinary DIS. The 
large momentum transfer means that the hard subpro- 
cess occurs on a spacetime scale much smaller than the 
bound state proton and is thus "embedded" in the pro- 
ton. The emerging hard-scattered partons, therefore, 
propagate through the proton's color field and may inter- 
act with it. Soft exchanges will dominate, due to the large 
coupling and the lack of suppression from hard gluon 
propagators. Therefore, the momenta of the hard par- 
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FIG. 1: 7*p — > Xp process with resummed gluon exchanges, 
and illustration of the factorization in Eq. (|4| of the amplitude 
into a hard and a soft part connected via an unintegrated 
gluon density function (UGDF). 



tons are essentially undisturbed, which is consistent with 
the fact that soft long distance interactions do not af- 
fect hard short distance ones. However, the exchange 
of color may change the color charges of the emerging 
partons such that the confining string-like field between 
them will have a different topology, resulting in a differ- 
ent distribution of the final state of hadrons produced 
from the string hadronization In particular, a region 
in rapidity without a string will result in an absence of 
hadrons there, i.e. a rapidity gap. This SCI model is 
very successful in describing data 6], but lacks a solid 
theoretical basis. 

Here, we present a new QCD-based model, which leads 
to effective color singlet exchange and thereby to diffrac- 
tive scattering. The model is inspired by the success of 
the SCI model, and may be seen as an explicit realization 
of the earlier attempt [7] to understand this soft gluon 
exchange in terms of QCD rescattering. 

As depicted in Fig. [T] the exchanged photon fluctu- 
ates into a quark-antiquark color dipole which interacts 
with the first gluon carrying a longitudinal momentum 
fraction xp. This defines the "hard" part of the process. 
Both the qq dipole and the proton remnant then have 
overall color octet charges and may interact through the 
exchange of a number of soft gluons with longitudinal 
momentum fractions x'j <C xp, '^x'^ = x' . This mul- 
tiple gluon exchange constitutes the "soft" part of the 
process and includes at least one exchange. Thus the 
overall exchange from both the hard and soft parts con- 
tains two gluons or more. The color factor of the total 
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soft exchange may combine to form an overall color sin- 
glet together with the first, hard gluon. In this case both 
the qq-system. and the proton remnant emerge as color 
singlets and hadronizc independently with a large sepa- 
ration in rapidity due to the dominance of small- a; gluons 
in the proton. Lacking a large momentum transfer to the 
proton, its remnant can recombine into a leading proton 
with close to the full beam momentum. The photon with 
a space-like virtuality has through the momentum ex- 
change xp been turned into a time-like hadronic system 
X of invariant mass Mx- 

In terms of the four-momenta q of the photon, P and 
P' of the initial and final proton, the important kincmat- 
ical variables are 
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where — —q^ and 



W' = {P + qf^^il-XB). (2) 
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In the forward limit of interest here, the momentum 
transfer t = {P' — P)^ along the proton line is small, 
|t| « A^. 

Let us now outline the model and the calculation of the 
diffractive structure function (see Q for details) . The in- 
volved momenta are specified in Fig. [T] We consider the 
asymmetric case where one hard gluon carries most of 
the longitudinal momentum transfer xp. Using cutting 
rules, we derive a factorization of the amplitude into a 
convolution of a hard part and a soft part. The hard 
part is treated in normal perturbative QCD. The soft 
part consists of any number of soft gluon exchanges, col- 
lectively in a color octet with total x' <^ xp. The soft 
exchanges are resummed, in the large Nc limit, to all or- 
ders in as- These soft gluons are not perturbative, since 
the strong coupling becomes large. We model them as 
interacting with quarks as perturbative gluons but with 
a non-perturbative coupling to be specified below. 

In the center-of-mass frame of the final state, i.e. the 
outgoing proton with momentum P', and the diffractive 
system X with momentum q' = ki + k2, where ki^j_ = 
— k2,_L = k^, the quark virtuality k'^ is the hard scale of 
the process, fij^ = k"^ = + k\ and is expressed in terms 
of its energy e and transverse momentum k± given by 



kl^z{l-z)Ml-ml, (3) 



where z is the fraction of the longitudinal momentum 
carried by the quark, and we consider the light quark 
mass limit uriq <C ■ 

In impact parameter space, with b conjugate to A_l, 
the total amplitude for the 7*p — ?> Xp process can be 
written as a convolution of the hard and soft subprocess 
amplitudes and a function V, which describes the distri- 
bution of gluons in the proton. 
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where S = = |A± + ^11- The function V is the 

Fourier transform of the unintegrated gluon distribution 
function (UGDF), and will be specified below. This fac- 
torization is schematically illustrated in Fig. [T] 

The amplitude for the hard subprocess 7*17 — >■ qq is 
decomposed into its longitudinal (L) and transverse (T) 
parts depending on the photon polarization. It includes 
the two possible couplings of the gluon to the qq pair, and 
can be Fourier transformed to impact parameter space 
with r conjugate to k'j_, the transverse momentum of a 
quark in the intermediate state. The hard amplitudes 
are then given by 
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where C — Sireq^/Tra^ / and Kq^i are Bessel func- 
tions. 

We now turn to the soft subprocess amplitude, which 
can be calculated order-by-order and then resummed. It 
is important to realize that these soft gluons carry non- 
perturbatively small momentum transfers. We deal with 
this by taking as(M) very small scales as a parameter, 
which we fix using the infrared-stable analytic pertur- 
bation theory (APT) [9!]. In the hmit that ^ 0, we 
use aff* = a^PT(Ai ^ Aqcd) ^ 0.7. The softness of 
the color-screening gluons with ^ xp implies that in- 
termediate particles are on-shell, and the dipole size r is 
frozen. Cutting the intermediate propagators, we pick up 
phase shifts originating from the hard amplitude, which 
depend on the soft momentum exchanges ^ . The dia- 
gram for one soft gluon exchange is a tree-level diagram, 
while two-gluon exchange leads to a loop integral. We 
calculate these contributions and perform the Fourier 
transforms with respect to A'^, where remarkably the 
second order diagram turns out to be the second term in 
a series that will exponentiate. This relies on the large 
Nc limit, where the color factors simplify to Cp — TpNc- 
We get 
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where A = 27ri CFas(A*soft)' A^Lft ^1 gluon 
virtuality, and we have defined 
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Summing over the number of soft gluons in the final state 
leads to exponentiation in impact parameter space, so 
that for the total soft subprocess amplitude we finally 
get 
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A similar expression was previously derived in the case 
of scalar Abelian gauge theory in Ref. Note, that 

M'*°^'(b,r) is independent of the photon polarization in 
the soft limit of small ^ . 

To describe the coupling of the gluons to the pro- 
ton, we use the framework of /cj^-factorization and gen- 
eralized (off-diagonal) UGDFs, which contain all infor- 
mation about the non-perturbative coupling of the glu- 
ons to the proton, and is based on a well-defined for- 
mal procedure for the transition from the parton level 
to the hadron level (see e.g. Ref. [HI). The cou- 
pling of a gluon to a quark is thus given by an off- 
diagonal UGDF Tf{xp,x',A]_,A'^_'^,n%), absorbing a 
factor CfCks (/isoft)/7r, and by convention a gluon prop- 
agator A]|^ into the UGDF in order to keep it regular 
as A5_ — >■ 0. The absorbed couphng asifJ-l^fi) corre- 
sponds to the coupling of a screening gluon with virtu- 
ality /igoft ^ A^ to a quark in the proton, whereas the 
coupling of the hard gluon to the qq dipole and to a quark 
in the proton is treated pcrturbatively at the hard scale 

Generalized parton distributions (GPDs) are not very 
constrained by data. We use a prescription for the gen- 
eralized UGDF, which was introduced in Ref. [T^], moti- 
vated by positivity constraints for GPDs [11] ■ This pre- 
scription works well in the description of recent GDF 
data on central exclusive charmonium production [l^ . 
and allows incorporating the dependence on the longitu- 
dinal momentum fraction and transverse momentum of 
the soft gluons in an explicitly symmetric way, 



on x' in the case x' ^ Xp: 
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which explicitly involves the soft x' dependence. Here 
J^g is the normal diagonal UGDF, which depends on the 
gluon virtuality, and which when integrated over this vir- 
tuality reduces to the well-known coUinear gluon PDF 
gix,^"^). The dependence of J-g on the virtuality is not 
theoretically well-known for small virtualities, and the 
UGDF is here modeled using the collinear gluon PDF 
together with a simple Gaussian Ansatz for the intrinsic 
transverse momentum dependence, 



^ xpg{xp, n%) x'g{x', nl^^^) fciAl), 
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where the factor y/xp is absorbed from the hard subpro- 
cess, and the Gaussian width po is the soft hadronic scale. 



corresponding to the transverse proton size 



1/Po. 



Note that this leads to an exponential i-dependence of 
the cross section ~ exp{Bp)t) with the diffractive slope 
Bd = 1/Po - 6-9 ± 0.2 GeV^ 

The second PDF in Eq. (fTTj) . associated with the soft 
gluon, is evaluated at very low scale and very small 
x' . For this PDF we can here introduce a function 
-Rg(a;', /Xg^fj) which is assumed to be slowly dependent 



xpJ^^ ~Rg{x\ )Jxpg{xp,p\)fG{A\ 
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The factor Rg, therefore, contains all the soft physics 
related with soft gluon couplings to the proton. It is in- 
terpreted as the square root of the gluon PDF at very 
small x' <^ Xp and some soft scale Atgoff This is a non- 
perturbative object, which contributes to the overall nor- 
malization and can be determined from data. The factor 
Rg in Eq. (|12p is analogous to the skewedness parame- 
ter Rg ~ 1.2 — 1.3, which accounts for the single logQ^ 
skewedness effect in off-diagonal UGDFs |16j| . As we will 
see below, the prescription (IT2|) is consistent with the 
HERA data for all available M|- and Q"^ . 

The model (fTTj) will lead to a linear dependence of 
the diffractive structure function on the gluon PDF, as 
compared to the quadratic dependence often encountered 
in two-gluon exchange calculations of DDIS p/Tj. This 
linear dependence is the same as in the SCI model, where 
a linear dependence describes both diffractive and non- 
diffractive events, and indicates a continuous transition 
between the two types of events. 

In terms of the UGDF in Eq. the factor V = 

V(b, r) of Eq. (g]) is given by 
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Straightforward calculations [3| lead to the final ex- 
pressions for the diffractive structure functions 

xpF^^^'^ ^SQ'^M\ dz{l~2z)z'^{l- zf\JL? 
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where S — J2q / {'^^'^ ^c) sums over light quark charges 
Cq, and 
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Let us briefly discuss the role of higher-order QCD cor- 
rections in the framework of our model. One should dis- 
tinguish corrections from the hard gluon emission due to 
radiation from the partons in the hard part, and correc- 
tions due to interactions between soft gluons in the soft 
part. 
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In the first case, additional s-channel gluons emit- 
ted from tfie hard scattering part can be described by 
DGLAP evolution. They play an important role for large 
invariant masses M|- ^ Q^, and will be considered be- 
low. 

In the second case, all interactions between soft gluons 
in the soft scattering part are absorbed into the UGDF 
and thus into the soft Rg factor, which enters the overall 
normalization. Only the number of soft gluon legs at- 
taching to the proton in the lower part of the diagram, 
and how they attach to the partons in the upper part 
are important. All long-distance interactions between 
the gluons are treated as part of the color background 
field in the proton, and do not affect the resummation 
procedure. 

One could also imagine interactions between the gluon 
from the hard part and one of the soft gluons. Such inter- 
actions contribute only in the symmetric case x' ^ xp, 
when all exchanged gluons are either soft or hard. The 
first case is unrealistic as it may happen only in the case 
of very small Mx and Q where QCD factorization does 
not apply. The second case is suppressed by a small a^. 
Such exchanges may be enhanced by large logarithms, 
leading to exchange of a BFKL pomeron at non-zero 
momentum transfer t. This is, however, a process with 
different kinematics and docs not contribute to forward 
diffraction. 

In the large- Mx, or /3 — > 0, limit the additional emis- 
sion of a gluon in the final state becomes important. This 
is dominated by the emission of a collinear gluon from the 
hard gluon, an emission which is enhanced by a large log- 
arithm. Such a gluon will be well-separated from the qq 
pair in momentum space, and will therefore contribute 
to building up a large Mx- We take this contribution 
into account through a gluon splitting, described using 
the DGLAP splitting function Pgg as (see e.g. [l^) 
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where tg is the gluon propagator and the integral is cut- 
off in the infrared by the effective gluon mass mg ~ 
Aqcd- The factor appears because the emitted 

gluon must contribute to the color singlet X system. 

The HERA data on the diffractive structure function 
[isj are given in terms of the reduced cross section. 



xpa. 
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where F^':^\xp, Q^, /3) is the diffractive structure func- 
tion integrated over t, the kinematical variable y = 
Q"^ I sxb < 1, and the center-of-mass energy of ep- 
cohisions at HERA is = 318 GeV. 

In Fig. [2] we show the comparison of the results of our 
model with the latest HERA data fl5| on the reduced 

cross section xpaP'''^\xp, /3,Q^) as a function of xp in 
bins of /3 and Q^. The figure shows our main result, 
calculated using the CTEQ6L1 gluon PDF [ll|, and also 



curves obtained using the older GRV94 PDF '20'l. The 
minimal factorization scale ^p is fixed to be IJ^^prnm — 
0.2 GeV^, which, together with Eq. ([3|) implies a minimal 
possible fraction of the quark longitudinal momentum 
Zinin in the integrals in Eqs. (|14I15I) . 

In these results, transverse polarization dominates in 
all bins. We find that the qq contribution alone is enough 
to describe all the data for /? > 0.2, below which the qqg 
contribution becomes significant. 

The fixed parameters in our model, which all take 
reasonable physical values, are an effective gluon mass, 
mg ~ Aqcd, used to regulate the infrared divergence 
in the qqg contribution, Eq. ([T7)) . and the soft coupling 
constant as(^gQfJ ~ 0.7. 

We also fit two physical quantities: the soft factor Rg, 
which absorbs the non-perturbative couplings of the soft 
gluons to the proton, and the constituent quark mass m^. 
For our model to be consistent, these two parameters 
should not depend strongly on the two large scales in the 
process, Mx and Q. 

We have found that Rg is close to unity for a wide 
range of scales. It does not depend at all on Mx, and 
only in the lowest Q bins is there a noticeable increase in 
the normalization Rg by at most a factor 4. Here, how- 
ever, QCD factorization becomes questionable and the 
conventional gluon PDFs are poorly known such that us- 
ing the GRV94 PDF instead of_CTEQ6Ll this increase is 
essentially removed ^] . Thus Rg is essentially an overall 
normalization factor close to unity, which contains un- 
known information on the density of soft gluons in the 
proton. 

The kinematics, e.g. Eq. ([3]), depends on an effective 
quark mass corresponding to a dynamic, dressed quark 
generated dominantly by softer gluon radiation which 
cannot be calculated theoretically. Therefore, the param- 
eter niq is fitted to data and found to have only a slow 
variation with the large scales Mx and Q; namely in the 
range 0.3-1.2 GeV consistent with mainly soft dynamics, 
as expected. 

The uncertainty in the PDFs is illustrated in Fig. [21 
where we show, in the bins of small and Mx, the 
results obtained using the GRV94 gluon PDF. In these 
bins GRV does better in reproducing the data, because 
of its larger gluon density at these small x and scales. 

In summary, in this Letter, we have presented a new 
QCD-bascd model of soft gluon exchanges in the fi- 
nal state, which describes data on the diffractive struc- 
ture function very well. The model is inspired by 
the phenomenologically successful Soft Color Interaction 
model ^] and on the work on such soft rescattering in 
DIS in Ref. [7|,[l3|- The full details of the model and the 
calculations are presented in Q. 

We have considered diffractive DIS, as this is where 
the most precise data are available, but the soft gluon 
exchanges arise due to the proton's color field and should 
thus be of a universal nature. Our model should therefore 
be applied to other processes, for example, diffraction in 
hadron-hadron collisions or diffractive vector meson pro- 
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FIG. 2: The reduced cross section xpa^^^\xp, P, Q^) as a function of xp for different values of Mx and . The latest 
ZEUS data , from diffractive deep inelastic scattering events with a large rapidity gap, compared with our model using the 
CTEQ6L1 (full line) and GRV94 (dotted line) parametrizations of the gluon density in the proton. 
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duction. It should also have effects on other observables, 
such as the underlying event at LHC, which is important 
to understand and describe. This is also borne out by 
the results \d\ from the SCI model, which has not only 
been able to describe all diffractive data from HERA and 



the Tevatron, but also other, non-diffractive data. 
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